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We present a measurement of the average value of a new observable at hadron colliders that is 
sensitive to QCD dynamics and to the strong coupling constant, while being only weakly sensitive 
to parton distribution functions. The observable measures the angular correlations of jets and 
is denned as the number of neighboring jets above a given transverse momentum threshold which 
accompany a given jet within a given distance AR in the plane of rapidity and azimuthal angle. The 
ensemble average over all jets in an inclusive jet sample is measured and the results are presented as 
a function of transverse momentum of the inclusive jets, in different regions of AR and for different 
transverse momentum requirements for the neighboring jets. The measurement is based on a data 
set corresponding to an integrated luminosity of 0.7 fb -1 collected with the DO detector at the 
Fermilab Tevatron Collider in pp collisions at y/s = 1.96 TeV. The results are well described by a 
perturbative QCD calculation in next-to-leading order in the strong coupling constant, corrected 
for non-perturbative effects. From these results, we extract the strong coupling and test the QCD 
predictions for its running over a range of momentum transfers of 50-400 GeV. 

PACS numbers: 13.87.Ce, 12.38.Qk 



Quantum Chromodynamics (QCD) predicts that the 
strong force between quarks and gluons becomes weaker 
when probed at high momentum transfers, correspond- 
ing to small distances. This property, referred to as 
asymptotic freedom, is derived from the renormalization 
group equation (RGE) [3-0]. The RGE does not pre- 
dict the value of the strong coupling a s , but it describes 
the dependence of a s on the renormalization scale /j,r, 
and therefore on the momentum transfer. Tests of per- 
turbative QCD (pQCD) and the property of asymptotic 
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freedom can be divided into tests of the validity of the 
RGE and determinations of the value of a s . By conven- 
tion, a s values extracted from data at different momen- 
tum transfers are evolved to the common scale (1r = Mz 
to allow comparisons between experiments. The current 
world average value is a s (M z ) = 0.1184 ± 0.0007 0]. 
The validity of the RGE is tested by studying the de- 
pendence of a s on the momentum transfer. At present, 
the RGE predictions have been tested in e + e~ anni- 
hilation, where a s results have been obtained for mo- 
mentum transfers up to 208 GeV Q ■ Attempts to ex- 
tract a s at higher momentum transfers have been car- 
ried out using inclusive jet cross section data in hadron- 
hadron collisions [!, 0] • These analyses methods require 
parton distribution functions (PDFs) of the proton at 
large scales as input. Since the main constraints on 
PDFs come from data at lower scales, the knowledge of 
PDFs at large scales is mainly based on the evolution 
according to the Dokshitzer-Gribov-Lipatov-Altarelli- 
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Parisi (DGLAP) evolution equations which use a s 
and the RGE as input. The a s results from inclusive 
jet cross section data at high momentum transfers can 
therefore not be regarded as tests of the RGE, since they 
are derived assuming its validity. 

In this Letter a new observable for hadron-hadron col- 
lisions is introduced and its average value is measured. It 
is related to the angular correlations of jets. In pQCD, 
this quantity is computed as a ratio of jet cross sections, 
which is proportional to a s . Since PDF dependencies 
largely cancel in the ratio, the extracted a s results are 
almost independent of initial assumptions on the RGE. 
Values of a s are extracted for momentum transfers be- 
tween 50 and 400 GeV. These provide the first test of the 
RGE at momentum transfers above 208 GeV. 

The analysis presented in this Letter studies the prop- 
erties of multi-jet production based on an inclusive jet 
sample in pp collisions at yfs = 1.96 TeV. While pQCD 
predictions for any cross section at a hadron collider de- 
pend on the PDFs, quantities with significantly reduced 
PDF sensitivity can be constructed. One class of such 
quantities is ratios of three-jet and dijet cross sections. 
Based on such ratios, one can exploit the high energy 
reach at hadron colliders to determine a s and to test the 
predictions of the RGE at previously unexplored momen- 
tum scales. A new observable is introduced, which probes 
the angular correlations of jets in the plane of rapidity 
y (loj and azimuthal angle <p. This observable measures 
the number of neighboring jets that accompany a given 
jet with transverse momentum (pt) with respect to the 
beam axis. The measured quantity Rar is the ensemble 
average over all jets in an inclusive jet sample of this ob- 
servable. The inclusive jet sample consists of all jets in a 
given data set, and these jets are hereafter referred to as 
"inclusive jets" . The measured quantity is given by 



y-^jct(pT) pr(i) (AR nbr 



(1) 



where N- ]ot (pT) is the number of inclusive jets in a given 
inclusive jet pr bin, and N^ T (AR, p T h ^ in ) is the number 
of neighboring jets with transverse momenta greater than 
-PTmim separated from the i-th inclusive jet by a distance 
AR within a specified interval AR m i n < AR < AR max 
with AR = ^/(Ay) 2 + (A<f>) 2 . For AR < tt, only topolo- 
gies with at least three jets contribute to the numerator 
of Eq. (P), in pQCD, and Rar is computed at lowest or- 
der as a ratio of three-jet (O(a^)) and inclusive jet cross 
sections (0(a 2 )). This ratio is proportional to a s . 

This measurement is based on a data set correspond- 
ing to an integrated luminosity of 0.7 fb -1 collected 
with the DO detector at the Fermilab Tevatron Collider. 
Rar{pt, AR, Pj^in) ^ s measured in an inclusive jet sam- 
ple at central rapidities \y\ < 1 for pt > 50 GeV, de 



measured triple differentially, as a function of inclusive 
jet pt, for different p^njm, an d m different AR regions. 
The Prmin requirements are 30, 50, 70, or 90 GeV, respec- 
tively, and the different AR intervals are 1.4 < AR < 1.8, 
1.8 < AR < 2.2, and 2.2 < AR < 2.6. For jets with 
-Rconc = 0.7, the lower limit of AR > 1.4 ensures that a 
jet does not overlap with its neighboring jets. The upper 
limit on AR is smaller than tt, so that contributing neigh- 
boring jets stem only from three- (or more) jet topolo- 
gies. The lowest Pj^in requirement is chosen to ensure 
that the jet energy calibration and the jet pt resolutions 
are well understood. The trigger efficiencies are high for 
jets with pt > 50 GeV in the inclusive jet sample. The 
requirement of \y\ < 1 implies that (\y\ + AR) < 3.6 over 
the whole analysis phase space. In this rapidity region 
jets are well-measured in the DO detector. The data are 
corrected for experimental effects and are presented at 
the "particle level," which includes all stable particles as 
defined in Ref. [ij. 

A detailed description of the DO detector can be found 
in Ref. [l3| ■ The event selection, jet reconstruction, and 
jet energy and momentum correction follow closely those 
used in recent DO measurements of inclusive jet, dijet and 
three-jet production rates [14l4l8l |. Jets are reconstructed 
in the finely segmented liquid- argon/uranium calorime- 
ter which covers most of the solid angle for polar angles 
of 1.7° <9< 178.3° [I3|. For this measurement, events 



are triggered by jet triggers. Trigger efficiencies are stud- 
ied as a function of jet pr by comparing the inclusive jet 
cross section in data sets obtained by triggers with differ- 
ent pt thresholds in regions where the trigger with lower 
threshold is fully efficient. The trigger with lowest pt 
threshold is shown to be fully efficient by studying an 
event sample obtained independently with a muon trig- 
ger. In each inclusive jet pt bin, events are taken from a 
single trigger which has an efficiency higher than 99%. 

The position of the pp interaction is determined from 
the tracks reconstructed using data from the silicon de- 
tector and scintilla ting fiber tracker located inside a 2 T 



solenoidal magnet [13j. The position is required to be 



fined by the Run II midpoint cone jet algorithm 11 1 



with a cone of radius R c 



0.7 in y and 4>. It is 



within 50 cm of the detector center in the coordinate 
along the beam axis, with at least three tracks pointing to 
it. These requirements discard (7-9)% of the events, de- 
pending on the trigger used. Contributions from cosmic 
ray events are suppressed by requiring the missing trans- 
verse momentum in an event to be less than 70% (50%) 
of the uncorrected leading jet pt if the latter is below 
(above) 100 GeV. The efficiency of this requirement for 
signal is found to be > 99.5% 14, 3|. Requirements on 
the characteristics of calorimeter shower shapes are used 
to suppress the remaining background due to electrons, 
photons, and detector noise that would otherwise mimic 
jets. The efficiency for the shower shape requirements 
is above 97.5%, and the fraction of background events is 
below 0.1% for all pt, as determined from distributions 
in signal and in background-enriched event samples. 



5 



10 — 



CL 

CC 



DC 



10 



10 r 




MSTW2008 PDFs 

NLOpQCD + 

non-perturb. correct. 



p nbr =30GeV 

Tmin 



D0 L = 0.7 fb" 




Vs = 1 .96 TeV 



p nbr =50GeV 

Tmin 




p nbr =70GeV 

Tmin 



* 2.2<AR<2.6 
■ 1.8<AR<2.2 

• 1.4 < AR < 1.8 




'p'- = 90GeV 

Tmin 



50 100 200 400 50 100 200 400 50 100 200 400 50 100 200 400 

P T (GeV) 



FIG. 1: (Color online.) The measurement of -Rah as a function of inclusive jet pr for three different intervals in AR and for 
four different requirements of p^mm- The inner uncertainty bars indicate the statistical uncertainties, and the total uncertainty 
bars display the quadratic sum of the statistical and systematic uncertainties. The theory predictions are shown with their 
uncertainties. 



The jet four-momenta reconstructed from calorimeter 
energy depositions are then corrected, on average, for the 
response of the calorimeter, the net energy flow through 
the jet cone, additional energy from previous beam cross- 
ings, and multiple pp interactions in the same event, but 
not for muons and neutrinos 3 EEi- The absolute 
energy calibration is determined from Z — > e + e~ events 
and the pt imbalance in 7 + jet events in the region 
\y\ < 0.4. The extension to larger rapidities is derived 
from dijet events using a similar data-driven method. In 
addition, corrections in the range (2-4)% are applied that 
take into account the difference in calorimeter response 
due to the difference in the fractional contributions of 
quark and gluon-initiated jets in the dijet and the 7 + 
jet event samples. These corrections are determined us- 
ing jets simulated with the PYTHIA event generator [2(| 
that have been passed through a GEANT-based detector 



simulation 21]. The total corrections of the jet four- 



momenta vary between 50% and 20% for jet pr between 
50 and 400 GeV. An additional correction is applied for 
systematic shifts in 



mm 



\y\ due to detector effects 
These corrections adjust the reconstructed jet energy to 
the energy of the stable particles that enter the calorime- 
ter except for muons and neutrinos. 

The differential distributions Rar(pt, AR,p T h ^ in ) are 
corrected for experimental effects. Particle-level events 
are generated with SHERPA [12] with MSTW2008LO 
PDFs J2p and with PYTHIA ^ with CTEQ6.6 
PDFs [24| and tune QW [25| . The jets from these events 
are processed by a fast simulation of the DO detector re- 
sponse. The simulation is based on parametrizations of 
jet pt resolutions and jet reconstruction efficiencies de- 



termined from data and of resolutions of the polar and 
azimuthal angles of jets, which are obtained from a de- 
tailed simulation of the detector using GEANT. 

The pr resolution for jets is about 15% at 40 GeV, 
decreasing to less than 10% at 400 GeV. To use the 
fast simulation to correct for experimental effects, the 
simulation must describe all relevant distributions, in- 
cluding the pt, y and AR distributions for the inclusive 
jets and the neighboring jets. The generated events are 
reweighted, based on the properties of the generated jets, 
to match these distributions in data. To minimize mi- 
grations between inclusive jet pt bins due to resolution 
effects, we use the simulation to obtain a rescaling func- 
tion in reconstructed pr that optimizes the correlation 
between the reconstructed and true values. The bin sizes 
in the pt distributions are chosen to be approximately 
twice the pt resolution. The bin purity after pt rescal- 
ing, defined as the fraction of all reconstructed events 
that were generated in the same bin, is above 50% for 
all bins. We then use the simulation to determine bin 
correction factors for experimental effects for all analy- 
sis bins. The correction factors are computed bin-by-bin 
as the ratio of Rar without and with simulation of the 
detector response. These also include corrections for the 
energies of unreconstructed muons and neutrinos inside 
the jets. The total correction factors for Rar using the 
reweighted pythia and SHERPA simulations agree typi- 
cally within 2%. The average factors, used to correct the 
data, are typically between 0.98 and 1.01, but never be- 
low 0.93 or above 1.03. The difference between the aver- 
age and the individual corrections is taken into account as 
an uncertainty which is split into two contributions. One 
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FIG. 2: (Color online.) The ratios of the Rar measurements and the theory predictions obtained for MSTW2008NLO PDFs 
and a s (Mz) = 0.118. The ratios are shown as a function of inclusive jet pr in different regions of AR (rows) and for different 
PTmin requirements (columns). The inner uncertainty bars indicate the statistical uncertainties, and the total uncertainty bars 
display the quadratic sum of the statistical and systematic uncertainties. The theory uncertainty is the quadratic sum of PDF 
and scale uncertainties. 



contribution corresponds to the systematic difference be- 
tween the two individual corrections, and the other one 
corresponds to the statistical fluctuations. The former 
is attributed to the model dependence and assumed to 
be correlated between the data points, while the latter is 
included in the statistical uncertainty of the results. 

In total, 69 independent sources of experimental sys- 
tematic uncertainties are identified, mostly related to jet 
energy calibration and jet pr resolution. The effects of 
each source are taken as fully correlated between all data 
points. The dominant uncertainties for the differential 
cross sections are due to the jet energy calibration (2— 
5)%, and the model dependence of the correction factors 
(2-3)%. Smaller contributions come from the jet pr res- 
olution (0.5-1.5)%, the jet ci resolution (0.5-2)%, and 
from the uncertainties in systematic shifts in y (0.5-1)%. 
All other sources are negligible. The total systematic 
uncertainties are between 2% and 6%. 

The results for Rar(Pt, Aii^Tmin) are displayed in 
Fig. [T] as a function of inclusive jet px, in different re- 
gions of AR and for different p^^ in - The values of pr at 
which the data points are presented correspond to the ge- 
ometric bin centers. A detailed documentation of the re- 
sults, including the individual uncertainty contributions, 
is provided in the supplemental material [26j . For a given 
AR region, and p^ in , Rar increases with pr up to a 
maximum value, above which it falls when approaching 



the kinematic limit. At fixed p?, Rar increases with AR 
and decreases with increasing Pxmm- At lower pr, Rar 
depends more strongly on Pxmin- For l ar g er PTmim b°th 
the px and the AR dependencies are stronger. 

The theory predictions for Rar, which are compared to 
the data, and which are later used to extract a s , are given 
by the product of the NLO pQCD results and correction 
factors for non-perturbative effects, including hadroniza- 
tion and underlying event. The non-perturbative correc- 
tions are determined using PYTHIA with tunes AMBT1 
and DW [27| . which use different parton shower and un- 
derlying event models. The hadronization correction is 
obtained from the ratio of Rar. on the parton-level (af- 
ter the parton shower) and the particle-level (including 
all stable particles), both without underlying event. The 
underlying event correction is computed from the ratio 
of Rar computed at the particle level with and without 
underlying event. The total corrections are defined as 
the combination of the corrections due to hadronization 
and the underlying event and they vary between +10% 
and -3% for tune AMBT1 and between -1% and -10% 
for tune DW. The results obtained with the two tunes 
agree typically within (2-4)% and always within 11%. 
The central results are taken to be the average values, 
and the uncertainty is taken to be half of the difference 
(given in Ref. (H). 

The NLO pQCD prediction is given by the ratio of an 
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inclusive three-jet cross section and the inclusive jet cross 
section both evaluated at their respective NLO. The nu- 
merator and the denominator both depend on the PDFs 
and most of the PDF dependencies cancel in the ratio. 
A residual PDF dependence remains, due to small differ- 
ences in the decomposition of the partonic subprocesses 
and a slightly different coverage of proton momentum 
fractions x in the numerator and the denominator. While 
the PDFs have no explicit a s dependence, their knowl- 
edge (i.e. PDF parametrizations) depends implicitly on 
a s due to assumptions on a s during the extraction pro- 
cedure. Therefore, the pQCD prediction for Rar has 
an explicit a s dependence stemming from the ratios of 
three-jet and inclusive jet matrix elements, and an im- 
plicit a s dependence due to the residual dependence on 
the PDFs. 

The NLO pQCD results are computed using 
FASTNL O I28l based on NLOJET++ [2^, [3y], in the MS 
scheme [31| for five active quark flavors. The calcula- 
tions use the next-to-leading logarithmic (two-loop) ap- 
proximation of the RGE and a s {Mz) = 0.118 in the 
matrix elements and the PDFs, which is close to the 
current world average value of 0.1184 [4]. The central 
choice fio for the renormalization and factorization scales 
is the inclusive jet px, Ur = (J-f = Mo = Pt, and the 
MSTW2008NLO PDFs [23| are used. 

The uncertainties of the pQCD calculations due to un- 
calculated higher order contributions are estimated from 
the fJ,R t F dependence. These are computed as the rela- 
tive changes of the results due to independent variations 
of both scales between (J,q/2 and 2/j,o, with the restriction 
of 0.5 < hr/ hf < 2.0. These variations affect the theory 
results by (3-9)%. The PDF uncertainties are computed 
using the up and down variations of the 20 orthogonal 
PDF uncertainty eigenvectors, corresponding to the 68% 
C.L., as provided by MSTW2008NLO. The i?AR_results 
obtained with the CT10 [H and NNPDFv2.1 [Hi PDF 
parametrizations agree with those for MSTW2008NLO 
typically within 1% and always within 3%. 

The theory results are compared to the data in Fig. [T] 
and the ratios of data and theory are displayed in Fig. [5] 
for all twelve kinematic regions in Ai? and Pj^n- The 
PDF uncertainties are (2-5)% and the scale uncertain- 
ties are typically (4-8)%. For higher p^mm = 70 > and 
90 GeV, the theoretical predictions are in good agreement 
with data and the ratios are independent of px, Ai?, and 
PTmin- Only for p^ in = 30 GeV, the predictions are sys- 
tematically below the data by (8-15)%. This might be 
caused by limitations of either the perturbative calcula- 
tion or the modeling of the non-perturbative effects at 

These Rar results are then used to determine a s and 
to test the two-loop RGE prediction for its running as a 
function of the scale px- In an initial study, the data are 
split into 12 subsets defined by the different (Ai?, Px£a n ) 
requirements. Assuming the RGE, the value of a s (Mz) 



TABLE I: The a s (Mz) results with their absolute uncertain- 
ties and the \ 2 values from the fits to the Rar data in each 



of the 12 kinematic regions, 
requirements. 



defined by the Pxmin and AR 



CTmiQ 


Ai? 


a J My) 


total 




y 2 /Nriof 

A, 1 ± 1 dot 










uncertainty 




30 GeV 


1.4- 


-1 8 


0.1290 


+0.0073 - 


0.0078 


6.9 / 11 


30 GeV 


1.8- 


-2.2 


0.1276 


+0.0078 - 


0.0049 


12.6 / 11 


30 GeV 


2.2- 


-2.6 


0.1249 


+0.0133 - 


0.0020 


15.3 / 11 


50 GeV 


1.4- 


-1.8 


0.1197 


+0.0089 - 


0.0061 


7.3 / 11 


50 GeV 


1.8- 


-2.2 


0.1168 


+0.0083 - 


0.0039 


14.1 / 11 


50 GeV 


2.2- 


-2.6 


0.1193 


+0.0076 - 


0.0043 


13.7 / 11 


70 GeV 


1.4- 


-1.8 


0.1168 


+0.0101 - 


0.0073 


4.9 / 9 


70 GeV 


1.8- 


-2.2 


0.1132 


+0.0069 - 


0.0047 


12.1 / 11 


70 GeV 


2.2- 


-2.6 


0.1156 


+0.0080 - 


0.0039 


16.8 / 11 


90 GeV 


1.4- 


-1.8 


0.1135 


+0.0084 - 


0.0087 


1.2 / 9 


90 GeV 


1.8- 


-2.2 


0.1136 


+0.0067 - 


0.0069 


9.7 / 9 


90 GeV 


2.2- 


-2.6 


0.1166 


+0.0099 - 


0.0083 


17.3 / 11 



is fitted to each of these subsets, and the corresponding 
X 2 values are determined that compare data and theory. 
Since each of these subsets covers a large inclusive jet px 
range, a violation of the RGE would be reflected in poor 
X 2 values. Furthermore, the comparison of the extracted 
a s (Mz) values allows to study the dependence of the re- 
sults on Ai? and/or p^ in . The data from kinematic 
regions m (Ai?, p^ b n r lin ) in which the a s (M z ) fit results 
are consistent with each other are then used in the subse- 
quent analysis. These data are split into 12 groups, each 
with the same inclusive jet px, combining data points 
for different (Ai?, Pj^ in ). For each group, a s is deter- 
mined at the corresponding px, and then evolved, using 
the RGE, to ^ R = M z . 

The a s extraction requires the theory predictions to be 
available as a continuous function of a s used in the ma- 



trix elements and PDFs. The global PDF fits |23j, |3J, |33 1 
do not provide the full a s dependence of their results, 
but only PDF sets at discrete values of a s (Mz), in in- 
crements of Aa s (Mz) = 0.001. A continuous a s (Mz) 
dependence for Rar is obtained, by cubic interpolation 
(linear extrapolation) of the theory results inside (out- 
side) the available a s (Mz) range. For the central results, 
we use MSTW2008NLO PDFs which cover the largest 
range of 0.110 < a s (M z ) < 0.130. The fits determine 
a s by using minuit [34J to minimize the % 2 function 35 1 
calculated from the differences between theory and data. 
All correlated systematic experimental and theoretical 
uncertainties are treated in the Hessian approach [35j |. 
except for the uncertainty due to the [Jlr^f dependence. 
The correlated statistical uncertainties are taken into ac- 
count via the covariance matrix. The a s results are ob- 
tained by minimizing % 2 with respect to a s and the nui- 
sance parameters for the correlated uncertainties. By 
scanning % 2 as a function of a s , the uncertainties are 
obtained from those a s values for which x 2 is increased 
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TABLE II: Central values and uncertainties due to different sources for the 12 cx s (pt) results obtained by combining the data 
at the same pr from all AR regions for PTmm = 50, 70, and 90GeV. All uncertainties are multiplied by a factor of 10 3 . 



Pt range 
(GeV) 


PT 

(GeV) 


OisijpT) 


50 - 


60 


55.0 


0.1353 


60 - 


70 


65.0 


0.1299 


70 - 


85 


77.5 


0.1232 


85 - 


100 


92.5 


0.1180 


100 - 


120 


110 


0.1154 


120 - 


140 


130 


0.1107 


140 - 


170 


155 


0.1070 


170 - 


200 


185 


0.1041 


200 - 


240 


220 


0.1050 


240 - 


280 


260 


0.1061 


280 - 


340 


310 


0.1049 


340 - 


450 


395 


0.0966 



total 
uncertainty 



statistical 



+ V.2 

-5.6 
+8.1 
-6.6 
+4.9 
-5.3 
+4.9 
-3.8 
+2.8 
-7.4 
+6.0 
-3.9 
+5.4 
-3.8 
+6.7 
-4.0 
+5.4 
-3.3 
+5.5 
-6.3 
+5.4 
-6.2 
+ 7.8 
-10.8 



experimental 


non-pert urb. 


Mb 1 WzOOoJN-LU 


nni? 




i . j 
correlated 


corrections 


uncertainty 


set 


variation 


+2.6 


+2". 5 


+ 1.3 


+0.2 


+5.4 


-2.8 


— 2.8 


— 1.2 


—0.4 


-0.8 


+2.3 


+2.1 


+ 1.2 


+0.3 


+6.1 


— 2 7 


— 2 4 


— 1 4 


— 1 3 


— 1 5 


+ 1.6 


+ 1.4 


+ 1.9 


+ 1.8 


+3.5 


— 3.2 


-1.0 


-1.0 


—0.9 


-3.9 


+2.8 


+ 1.0 


+2.1 


+ 1.0 


+3.0 


-2.4 


-2.2 


-1.1 


-0.0 


-1.4 


+2.1 


+0.3 


+ 1.0 


+0.0 


+ 1.4 


-2.4 


-0.4 


-5.0 


-3.7 


-3.1 


+2.8 


+0.4 


+ 1.5 


+2.0 


+4.7 


-2.2 


-0.4 


-2.5 


-0.0 


-1.9 


+ 1.6 


+0.1 


+0.9 


+ 1.5 


+4.9 


-3.0 


-0.3 


-0.8 


-0.0 


-2.2 


+2.5 


+0.7 


+0.3 


+3.0 


+ 5.4 


-2.1 


-0.4 


-1.5 


-0.0 


-2.9 


+2.5 


+ 0.6 


+ 1.0 


+0.8 


+4.5 


-2.3 


-0.3 


-0.2 


-0.6 


-2.3 


+ 1.0 


+ 1.0 


+0.3 


+0.0 


+ 5.3 


-3.2 


-0.8 


-0.7 


-3.3 


-4.2 


+ 1.6 


+0.4 


+0.3 


+0.6 


+ 5.0 


-2.3 


-0.3 


-0.6 


-3.3 


-4.3 


+ 1.9 


+0.1 


+0.2 


+0.0 


+ 5.3 


-5.9 


-1.0 


-0.9 


-3.3 


-4.7 



±2.8 
±4.2 
±0.6 
±0.8 
±0.6 
±0.6 
±0.5 
±0.5 
±0.3 
±0.6 
±1.0 
±5.4 



by one with respect to the minimum value. Fits, that 
determine a s (Mz) use the two- loop solution of the RGE 
to translate a s (Mz) values to the corresponding values 
of a s (pT) which enter the pQCD calculations for the dif- 
ferent pt bins. These a s (Mz) results are therefore as- 
suming the validity of the RGE. Those fits that extract 
(x s {pt) from a group of data points in the same px bin 
are almost independent of the RGE. A small dependence 
on the RGE enters only due to the residual dependence 
of the Rar predictions on the PDFs which use the RGE 
in their DGLAP evolution. Otherwise these o. s (pt) fit 
results are independent of the RGE. 

In the a s determination, we consider the correlations 
of the statistical uncertainties and all 69 sources of cor- 
related experimental systematic uncertainties. The the- 
ory uncertainties include the uncertainties of the non- 
perturbative corrections, the PDF uncertainties and the 
Hr,f dependence of the pQCD calculations. Following 
Refs. the uncertainty due to the /j,r,f depen- 

dence is computed by repeating the a s fit for different 
choices of p,R : F and the largest difference to the central 
result (obtained for /i^f = Pt) is taken to be the cor- 
responding uncertainty for a s . The a s fits are also re- 
peated for CT10 and NNPDFv2.1 PDFs, and the largest 
differences are quoted as "PDF set" uncertainty. The 
uncertainties from the scale variation and from the dif- 
ferent PDF sets are added in quadrature to the other 
uncertainties to obtain the total uncertainty. 

Before the central a s results are obtained, the con- 
sistency of the individual results for the 12 different 
(Ai?,p^ in ) regions, listed in Table UJ is tested. Assum- 
ing the RGE, the values of a s {Mz) are fitted to each 
of the 12 subsets, and listed in Table UJ together with 
the corresponding \ 2 values. All y 2 values are consis- 
tent with the expectations based on the number of de- 
grees of freedom (N do f), X 2 — N do f ± V^Ndof- This 
means that the RGE is consistent with the observed 
Pt dependence of a s (pT) over the studied pt range in 



all AR regions and for all p^mm- For the same p^ in , 
the a s (Mz) results for different AR regions are consis- 
tent with each other, i.e. there is no AR dependence. 
The a s (Mz) results are rather independent of Pj^; n for 
PTmin > 50 GeV. Only the a s (M z ) results for the lowest 
30 GeV, are significantly higher. 



requirement, p^ in 
As mentioned earlier, at lowest 



nbr 



Pf 



limitations of the 



perturbative calculations or the non-perturbative models 
may become visible. The data with p^mm = 30 GeV are 
therefore excluded when the final results of this analysis 
are determined. 

All remaining data points with the same pt (from all 
three AR regions and for p T h ^ in = 50, 70, and 90 GeV) 
are combined to fit ol s {pt), at the pt value correspond- 
ing to the geometric center of the bin. This is done for 
all 12 different pt bins in the range 50 < pt < 450 GeV 
and the results are listed in Table [IT] and displayed in 
Fig. [3] (a). Using the RGE, the individual results are 
then evolved to p,R = Mz, and shown in Fig. [3] (b). 
These a s results from Rar, extracted using NLO pQCD, 
are in good agreement with our previous results from in- 
clusive jet cross section data [36j, extracted using NLO 
plus 2-loop contributions from threshold corrections [3i| , 
and with the results from a reanalysis of event shape 
data from the ALEPH experiment at the LEP e + e~ col- 
lider, extracted using NNLO calculations (53]. A com- 
bined fit, using the same data set integrated over pt, 
and for MSTW2008NLO PDFs, gives the a s (M z ) result 
listed in Table El] The results obtained for CT10 PDFs 
(a s {M z ) = 0.1189) and NNPDFv2.1 {a s (M z ) = 0.1167) 
are used to define the uncertainty due to the PDF set. 
This result is in good agreement with our previous result 
of a s (Mz) — 0.116ltaoo48' obtained from inclusive jet 
cross section data at pt < 145 GeV and the world 
average value The RGE prediction for this result is 
displayed in Fig. [3] (a). The new « s (pt) results from Rar 
are well described by the RGE prediction including the 
region 208 < p,R < 400 GeV, in which the RGE is tested 
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FIG. 3: (Color online.) The strong coupling a s at large mo- 
mentum transfers, Q, presented as a„(Q) (a) and evolved to 
Mz using the RGE (b). The uncertainty bars indicate the 
total uncertainty, including the experimental and theoretical 
contributions. The new a s results from -Rah are compared 
to previous results obtained from inclusive jet cross section 
data [3(| and from event shape data [4(J. The a s (Mz) result 
from the combined fit to all selected data points (b) and the 
corresponding RGE prediction (a) are also shown. 



for the first time. 

In summary, a measurement has been presented of 
a new quantity Rar which probes the angular correla- 
tions of jets. Rar is measured as a function of inclusive 
jet pt in different annular regions of AR between a jet 
and its neighboring jets and for different requirements on 
the minimal transverse momentum of the neighboring jet 
PTmin- The data for pt > 50 GeV are well-described by 
pQCD calculations in NLO in a s with non-perturbative 
corrections applied. Results for a s (pT) are extracted us- 
ing the data with Pxroiu — 50 GeV, integrated over AR. 
The extracted a s {pr) results from Rar are, to good ap- 
proximation, independent of the PDFs and thus inde- 
pendent of assumptions on the RGE. Therefore, these 
a s results are the first to provide a test of the RGE at 
momentum transfers beyond 208 GeV. The results are in 
good agreement with previous and consistent with the 
RGE predictions for the running of a s for momentum 
transfers up to 400 GeV. The combined a s (Mz) result, 
obtained using the data with p^mm — 50 GeV (integrated 
over AR and p T ), is a s {M z ) = 0.1191^;^, in good 



agreement with the world average value |4 
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